Decay of highly-correlated spin states in a dipolar-coupled solid 
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We have measured the decay of NMR multiple quantum coherence intensities both under the 
internal dipolar Hamiltonian as well as when this interaction is effectively averaged to zero, in the 
cubic calcium fluoride (CaF2) spin system and the pseudo one-dimensional system of fluoroapatite. 
In calcium fluoride the decay rates depend both on the number of correlated spins in the cluster, 
as well as on the coherence number. For smaller clusters, the decays depend strongly on coherence 
number, but this dependence weakens as the size of the cluster increases. The same scaling was 
observed when the coherence distribution was measured in both the usual Zeeman or z basis and the 
x basis. The coherence decay in the one dimensional fluoroapatite system did not change significantly 
as a function of the multiple quantum growth time, in contrast to the calcium fluoride case. While 
the growth of coherence orders is severely restricted in this case, the number of correlated spins 
should continue to grow, albeit more slowly. All coherence intensities were observed to decay as 
Gaussian functions in time. In all cases the standard deviation of the observed decay appeared to 
scale linearly with coherence number. 
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I. INTRODUCTION 

While there have been several proposals put forward 
for scalable quantum computing architectures, experi- 
mental realizations have been limited to a handful of 
qubits at most. Maintaining coherence as the size of the 
system Hilbcrt space increases remains extremely chal- 
lenging. It is essential to understand how decoherence 
rates in different physical systems scale as a function of 
system size. There have been a number of theoretical 
investigations on the scaling behaviour of decoherence 
0- 13- H- These general models have typically been 
based on the spin-boson model of Leggett and co-workers 
P| . System-specific scaling laws have also been proposed 
for a few physical implementations (e.g. 0,0)- 

Palma et al. Q showed that for a multi-qubit quantum 
register, the decay of particular off-diagonal elements of 
the system density matrix depended on the Hamming 
distance / between the two states. In the case of inde- 
pendent, uncorrelated noise, the decay was of the form 
exp(— /r(t)) while for the case of correlated noise the 
decay was exp(— f 2 T(t)) where T(t) corresponds to the 
single qubit decay. 

Suter and coworkers recently published an experimen- 
tal study of the decay of multi-spin states using nu- 
clear magnetic resonance (NMR) [§[ . They used multiple- 
quantum (MQ) NMR experiments to create correlated 
multi-spin states in a powdered sample of the plastic 
crystal adamantane, and observed the rate at which these 
states decay during evolution under the internal dipolar 
coupling of the spins. They observed that the decay rate 



increased as a square root of the estimated number of 
correlated spins in the cluster. 

A theoretical analysis of their experimental results has 
been published recently by Fedorov and Fedichkin 
Neglecting the flip-flop (XY) terms of the dipolar inter- 
action, they obtained the following expression for the de- 
cay of multiple quantum coherence states, in the short 
time, large spin limit 



S n {t) = pcxp (— an 2 t 2 ) + (1 — p)exp ^— ~^ a tj 



(1) 
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where p = l/N(J2j djk) 2 / d 2 k is a correlation param- 
eter, a is proportional to the second moment of the line- 
shape, n is the coherence number and N is the number 
of correlated spins. The first term depends strongly on 
the coherence number and indicates that the spins in 
the multi-quantum state experience a correlated mean 
field. The second term does not depend on the coher- 
ence number, but only depends on the number of spins 
in the cluster, and indicates that the fields experienced 
by the different spins are uncorrelated from each other. 

These equations agree well with the measured data at 
short times, but deviate at longer times. The correlation 
parameter p appeared to remain constant with increasing 
spin system size. However since adamantane is a plastic 
crystal in which the molecules are undergoing rapid rota- 
tional motions, the intramolecular dipolar couplings are 
averaged to zero, and the only residual couplings are mo 
tionally averaged intermolecular dipolar couplings. This 
motional averaging modulates the contributions in Equa- 
tion 1 , making the experimental results obtained peculiar 
to this class of samples. 

In this paper we expand on these preliminary stud- 
ies. Our test systems are the cubic lattice of 100% abun- 
dant 19 F spins in a single crystal of calcium fluoride, and 
the pseudo one dimensional spin chains of fluoroapatite 
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(FAp). These are both rigid crystals, and molecular mo- 
tions will not affect the results obtained. In addition to 
characterizing the decay of the system under the dipo- 
lar Hamiltonian, we measure the decay rates obtained 
when we suspend evolution under the dipolar Hamilto- 
nian. We also repeated the experiments, encoding the 
multi-spin coherences in a different basis, and observed 
the resulting scaling behaviour. 

The standard MQ NMR experiment is shown in Figure 
1 0,^3- Assuming that the system is closed, the final 
density matrix is given by 



Ul Q U ev U ( (4>)U DQPi Ujjq U\ (</>) U\ v U DQ (2) 



4> is the phase angle and C, — {x, z} is the rotation 
axis for the multiple-quantum encoding. Collective ro- 
tations about each axis yield coherence order distri- 
butions in the corresponding basis. The propagator 
Udq = cxp (— UHdq) represents evolution under the 
double quantum (DQ) Hamiltonian given by 



is 

j<k 



D jk {erf at 



°3 a k 



} 



(3) 



where the dipolar coupling constant Djk between spins j 
and k is 



7 



(1 -3 cos 2 8 



(4) 



7 is the gyromagnetic ratio, rjk is the distance between 
spins j and fc, and 6jk is the angle between the external 
magnetic field and inter-nuclear vector r^. This effec- 
tive DQ Hamiltonian is created by multiple-pulse NMR 
techniques that toggle the dipolar Hamiltonian to create 
the appropriate zeroth order Average Hamiltonian. In 
Equation El above we have assumed that the experimen- 
tal implementation of —TLdq is perfect. The encoding of 
the coherence orders is performed by the collective rota- 
tions Uc_(4>) = cxp (^—i4>J2j a (j) ■ The observed signal is 
given by the overlap 



S(t) = Tr 
= Tr 



(5) 
(6) 



DQ 



where pf QC {(j)) = U ev lI (: ((f))p lD QU^((j))U^ v , and p 

UdqPiUdq, and the evolution U ev is defined below. 

Following this evolution the density operator of the 
spin system in the Zeeman basis contains off-diagonal 
terms of the form 



■a 1 a 2 



_n+l 



■ (T 



(7) 



We are interested in the properties of such coherences 
in the system. In the experiments here, we create 
a distribution of states with different coherence orders 



FIG. 1: A generic MQ NMR experiment. A tt/2 pulse is 
applied at the end of the pulse sequence for signal acquisition. 



M = (m — n) , and spin numbers K = n + m + p. After 
this step, we either allow the spin system to evolve under 
the internal dipolar interaction or suspend evolution of 
the spin system by applying a time suspension sequence. 
In the first case, U ev — exp (—iHot), where the dipolar 
Hamiltonian is 



Ti-D = D ]k \<y ]z a kz - -(<Tj+(Jk- + &j-&k+) \ ■ 
j<k L J 



(8) 



In the second case U ev = X, if the time suspension se- 
quence is perfect and the system is completely isolated. 
Decays observed during this time could be the result of 
errors in the control, or couplings to an environment. 

The experiments were performed at room temperature 
at 2.35 T (94.2 MHz, 19 F), using a Bruker Avance spec- 
trometer and home built probe. The samples used were 
a 1 mm 3 single crystal of CaF2 with Ti ~ 7s, and a 
crystal of fluorapatite (FAp) with T\ ~ 200 ms. The 
FAp crystal is a mineral crystal specimen from Durango, 
Mexico. High power 0.5 p,s ir/2 pulses were used. The 

rotation U z {cf>) — exp ( — icf> cr Z j\ was used to encode 
coherence number in the Zeeman or z basis while the 
rotation U x ((f>) = exp (— itfiHj ° "xj ) was used to encode 

coherence orders in the x-basis Il3l [13. The phase (0) 
was incremented from to Air with A(f> = ^ to encode 
up to 32 quantum coherences for every experiment. A 
fixed-time point corresponding to the maximum inten- 
sity signal was sampled for each <p value, and then was 
Fourier transformed with respect to <p to obtain the co- 
herence order distribution, as seen in Figure 

As the evolution time r under the DQ Hamiltonian 
increases (also referred to as the MQ growth time), pro- 
gressively more spins are correlated into the MQ states. 
Table 1 shows the size of the spin system for each MQ 
growth time, as estimated by the method of Baum et 
al. 0, and used in 0. A log-log fit indicates that the 
number of spins is increasing as N(= n + m + q) ~ t 2 . 



MQ growth time 


System size (N) 


43.3 p,s 


6 


86.8 ijs 


8 


130.2 fis 


12 


173.6 /is 


28 


217 {is 


36 


260.4 lis 


66 


303.8 lis 


96 



Table 1. Effective size of the spin cluster as a function of the 
MQ growth time, using the model in jlfl |. 
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FIG. 2: Zeeman basis coherence order distribution in CaF2 
with t=0, for an MQ growth time r=303.8 /is 
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FIG. 3: Decay of intensity for the different coherence orders 
for MQ growth times r = 130.3^is, and r = 303.8/is in CaF2. 
The solid lines are Gaussian fits to the data. 



II. EVOLUTION UNDER THE SECULAR 
DIPOLAR HAMILTONIAN 

The MQ states are not stationary under the internal 
dipolar Hamiltonian of the spins, and evolve as a function 
of the dipolar interaction time. This uncompensated evo- 
lution leads to imperfect refocusing during —Hdq- Fig- 
ure OH shows the intensities of various coherence orders in 
the z basis as a function of dipolar evolution time for two 
different MQ growth times. The signal intensity appears 
to decay as a Gaussian function in time. We have fit the 
data to exp(— t 2 /2Tj) to extract the effective decay times 
(Td = standard deviation) for each coherence order. The 
Gaussian shape of the decays indicates that the underly- 
ing process is consistent with a time-invariant dispersion 
of fields (due to the spins) having a normal distribution 

m 

Figure 01 shows the decay times Td of the different z 
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FIG. 4: Effective decay times under the dipolar Hamiltonian 
of the various coherence orders for different MQ growth times 
inCaF 2 . 



basis coherence orders as a function of the MQ growth 
times. Four features are evident in the data: (i) the Td 
are seen to depend linearly on the coherence number; 
(ii) for short r, the Td are seen to depend strongly on 
the coherence number; (iii) this dependence on coherence 
number weakens significantly with increasing t and (iv) 
the incremental change in Td with increasing r decreases. 

We performed a linear fit of Td versus coherence num- 
ber for each of the MQ growth times as 



T d = T d {0) -K-n 



(9) 



where Td(0) is the decay of the zero-quantum intercept, 
and k is the slope. Figure shows the dependence of 
ln(Td(Q)) and ln(— k) on ln(iV) the size of the spin system 
from Table 1, as well as the best linear fits obtained. 
We get ln(T d (0)) = 3.94 - 0.57 • ln(N) and Ixi(-k) = 
2.75- 1.21 -hx(N). Thusr d (0) w and k k -B/N, 

where A = 51.4 and B = 15.6. We can therefore express 
the scaling behavior of Td as 



T d 



A 
N 



B ■ n 

N 



(10) 



We have repeated the experiment encoding the coher- 
ence orders in the x-basis instead, and obtained identical 
scaling behaviour. The decays were once again observed 
to be Gaussian. In Figure [fj] we show the results of a 
multi-dimensional experiment in which we correlate the 
z and x basis decay times for an MQ growth time r = 
217.2 /is. It should be noted that the dipolar Hamil- 
tonian in a strong external magnetic field is anisotropic 
(see Equation[SJ), but this is not reflected in the observed 
decay times. 



4 











~ • - ~ 








ln(T d {0)) 
























■ 




■ " •» ^ 








"* " - - ■ In(-K) 






■ ■- „, 















2.5 3 3.5 4 4.5 

ln(size of spin system) 



FIG. 5: Dependence of ln(T d (0)) and 1h(-k) on ln(jV) the 
size the spin system, as well as the best linear fits obtained. 
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FIG. 7: Effective decay times under the 48-pulse sequence of 
the various coherence orders for different MQ growth times in 
CaF 2 . 



T=217.2iiS 




coherence number 



FIG. 6: Effective decay times under the dipolar Hamiltonian 
for correlated x and z bases coherence orders in CaF2. 
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FIG. 8: Ratio of decay time under the 48-pulse sequence to 
decay time under the dipolar evolution in CaF2. 



III. EVOLUTION UNDER A TIME 
SUSPENSION SEQUENCE 

We then attempted to suppress evolution of the dipolar 
Hamiltonian using a time-suspension pulse sequence that 
implements (approximately) the identity operator on the 
spin system [15J. In the ideal experiment, we should see 
no decay of the spin coherences due to dipolar couplings 
within the spin system. The cycle time of the 48-pulse 
time suspension sequence used here was 132.48 [is. The 
change in intensity of the z-basis coherence order was 
measured as a function of the number of loops of the 
48-pulse sequence. 

Once again we observed Gaussian decays as a function 
of time for the coherence intensities, indicating that the 
underlying noise has a long correlation time. We fit the 
data to a Gaussian and extracted the Td decay times. 
Figure [7] shows the effective T,j's of the z-basis coherence 



orders, under the time-suspension sequence. We see that 
the data is qualitatively identical to that obtained in the 
case of dipolar evolution, showing the same features dis- 
cussed above. This correspondence was repeated in the 
x-basis data, as well as in the correlated z and x basis 
encoding experiments. 

Figure [S] shows the ratio of the z-basis multiple quan- 
tum 2Vs during evolution under the 48-pulse time sus- 
pension sequence to the measured during dipolar evo- 
lution. The ratio was measured to be around 70, and 
appeared to be independent of the size of the spin corre- 
lations involved. Given the uniform scaling obtained, we 
did not repeat the linear fits. However it is clear from 
the uniform scaling that both A and B are scaled by the 
same factor of around 70. 
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FIG. 9: Effective decay times under the dipolar Hamiltonian 
in ffuoroapatite for MQ growth times r = 480/is and r = 
1200/is 



IV. QUASI-1D FLUOROAPATITE SYSTEM 

We also measured the decay of multiple quantum co- 
herences under the dipolar Hamiltonian in fluoroapatite. 
Floroapatite is a quasi-one dimensional spin chain, as 
the distance between spin chains is about three times 
larger than the distance between adjacent spins in the 
chain |16|. The one dimensional spin chain with nearest 
neighbor double quantum Hamiltonian is exactly solv- 
able [rH ITsI ] , and it has been shown that starting from a 
thermal equilibrium state, only zero and double quantum 
coherences are produced, even as higher order multi-spin 
states are created. The presence of higher order coher- 
ences indicates the importance of next-nearest neighbor 
and other distant couplings. Consequently it has been 
observed that higher order coherences grow very slowly 
in this system |19j. 

Consistent with the one dimensional nature of the spin 
system, the growth of the spin clusters is much slower in 
this case. More importantly, we see that the character of 
the decays does not change as a function of the the MQ 
growth time. A strong dependence on coherence number 
is observed for both short and long MQ growth times. 
However, even at long times, the number of correlated 
spins is still small in this system. 



V. DISCUSSION 

It is important to ensure that the observed decays are 
not an artifact of the encoding used to make the mea- 
surement. For example, we need to examine if the phases 
introduced into the system by the encoding U^((f>) lead 
to imperfect refocusing in our experiments. We mea- 
sured the decay time for the particular situation where 
U((<fr) = U^(ip) = I, as a function of the MQ growth 



time. The observed scaling for the case of dipolar evolu- 
tion was T d w 22.7/VN and T d ps 2276/ViV for the 48- 
pulse time suspension sequence. These times are slightly 
shorter than the corresponding values of ld(0) for the 
zero quantum coherence. This is reasonable as we are 
measuring the aggregate decay of all the MQ states here, 
and suggests that the encoding step is not responsible for 
the decay rates observed. 

The similarity between the dipolar evolution and the 
time suspension data suggests that the dominant source 
of noise in the time suspension data are residual dipo- 
lar coupling terms that are not effectively averaged out, 
as these would be expected to scale identically in the 
two cases. However, it is worth examining an alternative 
model in which the decay is due to coupling to an exter- 
nal environment. Fedorov et al. have also calculated the 
effect of a large, bath coupled to the multiple quantum 
states, and get a similar solution to that obtained under 
an inter-spin Ising coupling. Thus both theories - resid- 
ual dipolar errors and the presence of external spins yield 
identical scaling behaviours and it is not possible a priori 
to distinguish between theses two models on the basis of 
the data here. However we can examine the physics of 
the system under study to understand the origin of the 
decays. 

In principle the environment that the 19 F nuclear spins 
are coupled to could be lattice phonons or other spins — 
both electron and nuclear — that are present in the sys- 
tem. In order to effectively relax the nuclear spins, the 
phonons would need to be resonant with the Larmor fre- 
quency of the spins (94.2 MHz in these experiments). 
The Debye Temperature of calcium fluoride is 510 K, so 
it can be assumed that the spins are fixed in a rigid lat- 
tice in these room temperature measurements, and that 
phonons do not play a significant role in the relaxation 
of the spins. In addition to sparse paramagnetic impu- 
rities (responsible for T\ relaxation), we know that the 
CaF2 system contains 43 Ca which has spin 7/2 and is 
0.13% abundant. There could possibly be other spin de- 
fects in the crystal such as protons — albeit at much lower 
concentrations. 



Paramagnetic impurities 

We also measured the decay rates under the dipolar- 
evolution in a second crystal with a longer T± relaxation 
time pa 110 s, corresponding to a smaller concentration 
of paramagnetic impurities, and found no difference in 
the experimentally observed TVs. In Mn-doped CaF2 
with a T\ of 700 ms, Tse and Lowe estimated that the 
concentration of impurities (N e ) was approximately 5.6 x 
10 17 cm~ 3 [2Cj |. yielding an average impurity separation 



of 15 nm. Since \jT\ cx N e , a 7 s T\ corresponds to 
an impurity concentration of 5.6 x 10 16 cm -3 and an 
average separation of 32.4 nm, while for T\ = 110s, we 
get a concentration of 3.56 x 10 15 cm -3 and a separation 
of 81.3 nm. 



6 



The magnetic moment of a S=5/2 paramagnetic impu- 
rity is about 3500 times larger than that of the fluorine 
nucleus. The dipolar coupling between such an impurity 
and a 19 F nucleus becomes comparable to the strongest 
19 F- 19 F coupling (a = 2.73 A) at a distance of about 4 
nm. This corresponds to an interaction strength of about 
15 kHz, and dipolar correlation time of about 66 /is. This 
is approximately the size of the frozen core of fluorine 
spins around each impurity site, in which the "flip-flop" 
or XY terms of the dipolar Hamiltonian are suppressed. 
The electron-nuclear interaction further drops to 10 % of 
the inter-nuclear coupling at a distance of 8.9 nm. Thus, 
assuming Mn impurities, for a T\ = 7 s crystal, 84 % of 
the nuclear spins experience a hyperfine field that is less 
than 1.5 kHz. The correlation time of the hyperfine field 
seen by the bulk of the nuclear spins is much too long to 
explain the observed decays. 

43 Ca spins 

The magnetic moment of the calcium spins is about 
half that of the fluorine spins, and the Ca-F spacing is 
about 2.36 A, just a little shorter than the F-F distance. 
Thus we expect the strongest Ca-F coupling to be 11.6 
kHz. However, given the low natural abundance of the Ca 
spins (0.13 %), very few nuclear spins see a Ca coupling 
of this strength. Another possibility is that the Ca-Ca 
dipolar coupling leads to mutual spin flips, that reduces 
the efficiency of the 48-pulse time suspension sequence. 
The mean separation between Ca spins is 4.1 nm, and 
the average 43 Ca- 43 Ca coupling is about 2.2 Hz. Since 
the Ca-F coupling is decoupled on the time-scale of one 
cycle of the 48-pulse sequence, which is about 44 /is, the 
2 Hz Ca-Ca coupling is too weak to affect the efficiency 
of the decoupling sequence. Thus we do not expect the 
43 Ca spins to be the source of the observed decays. 

Errors in control 

We have assumed to this point that the implementa- 
tion of all the pulse-sequences has been perfect. The 
propagators Udq and U D q are idealized propagators 
for the NMR pulse sequences used, and typically cor- 



respond to the zeroth order Average Hamiltonian of the 
sequence. In reality the presence of higher order terms in 
the Magnus expansion result in U^U^ ± 1. More 
importantly, errors in the implementation of the two 
sequences — the presence of phase transients during the 
leading and falling edges of the pulses, as well as errors 
in the setting of the ir/2 pulse lead to imperfections in the 
refocusing. This imperfect refocusing is however unlikely 
to be the source of the observed decay in either experi- 
ment. In the dipolar evolution experiment, the strength 
of the dipolar coupling is much stronger than any residual 
error terms in the propagator, and it is these couplings 
that determine the decay rate. 

In the time-suspension experiments, the 48-pulse se- 
quence averages the dipolar interaction to zero to second 
order in the Magnus expansion. Assuming perfect im- 
plementation, the leading error terms are likely to be 
the second-order dipolar-offset term and the second or- 
der offset term. Our results indicate that the source of 
the observed signal decay in the time suspension exper- 
iments is the residual errors in the zeroth-order Aver- 
age Hamiltonian, which are proportional to TLd, rather 
than second-order terms like [Hd^i), TLd^)}, or other 
higher-order terms. The configurational space accessible 
by a single spin-flip is much smaller than that accessible 
by two spin-flips, especially at larger spin numbers. This 
suggests that we should expect significantly different scal- 
ing of decay rates for the two types of processes. Given 
the identical scaling behaviour observed in the two cases, 
we can conclude that the the dominant error in the time- 
suspension sequences has the form of TLd (possibly tog- 
gled about some arbitrary axis). The more likely source 
of error in the control is imperfect implementation of the 
sequence. Phase transients during the pulses, or small er- 
rors in the setting of the 7r/2 pulse can accumulate over 
the hundreds of pulses that are applied in this experi- 
ment, and the residual error terms will look like mod- 
ulated dipolar interactions, and consequently will scale 
the same way, though with a reduced strength [2lJ . 
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